Magnetic structure and orbital ordering in BaCoC>3 from first-principles calculations 
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Ab initio calculations using the APW+lo method as implemented in the WIEN2k code have been 
used to describe the electronic structure of the quasi-one-dimensional system BaCoOa. Both, GGA 
and LDA+U approximations were employed to study different orbital and magnetic orderings. GGA 
predicts a metallic ground state whereas LDA+U calculations yield an insulating and ferromagnetic 
ground state (in a low-spin state) with an alternating orbital ordering along the Co-Co chains, 
consistent with the available experimental data. 

PACS numbers: 71.15.Ap, 71.15.Mb, 71.20.-b 
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I. INTRODUCTION 

In the past few years, the main interest in Co oxides has 
been their applications in solid state fuel cells. 1 Further- 
more, some mixed valence Co compounds showed inter- 
esting colossal magnetoresistance (CMR) effects^ Very 
recently, a renewed interest in Co oxides started after the 
discovery of the first Co-based superconductori 3 . Due to 
these two phenomena, CMR and superconductivity, Co 
oxides are of great interest from a physical point of view. 
Most of the interesting physical properties in transition 
metal oxides lie in the interplay between electronic struc- 
ture and magnetic properties. New experimental tech- 
niques allow to investigate the coupling between mag- 
netic and orbital degrees of freedom, which plays a major 
role in this type of materialsA^ 

BaCoOa is a transition metal oxide with interesting 
properties, which have been studied in recent yearsA&£ 
In its highly anisotropic structure, face-sharing Co06 
octahedra form chains that are likely to produce one- 
dimensional effects in its magnetic properties. The 
hexagonal structure, however, may lead to magnetic frus- 
tration when the Co atoms perpendicular to the chains 
(on a triangular arrangement) should have an antiferro- 
magnetic coupling. 

The experimental studies on BaCoOa have established 
that the Co ions are in a low-spin state (S=l/2)jl as 
to be expected for a Co 4+ ion, with the configuration 
t2g e< g m the case of octahedral symmetry. In the temper- 
ature range from 70 to 300 K, the material is found to 
be semiconducting, with conduction occuring through n- 
type carriers^ The origin of the gap is not yet resolved. 8 
Possible reasons are a Mott-Hubbard-like transition, an 
Anderson localization due to some structural disorder or 
a "Peierls"-dimerization of the Co chains, but none of 
them has been confirmed experimentally. 

First measurements reported a Neel temperature 
Tjv=8 K^A but recent experiments reveal a more com- 
plex magnetic structure^ beyond the simple antiferro- 
magnetic behavior. Paramagnetism is found above 250 
K, where a coupling constant of J/ks = 10 K is estimated 



by fitting to a ID Heisenberg model. For temperatures 
between 70 and 250 K, ferromagnetic coupling is pre- 
dominant, whereas for T<70 K antiferromagnetic cou- 
plings become more important. Unfortunately no neu- 
tron diffraction data are available and the magnetic or- 
der at very low temperatures is not known. Also the rel- 
ative strengths of intra-chain and inter-chain magnetic 
couplings are still uncertain. 

In this paper, we try to shed some light into the 
electronic structure of BaCoOa by means of ab initio 
full-potential augmented plane wave plus local orbitals 
(APW+lo) calculations^ carried out with the WIEN2k 
codeiiSii 3 . We study some possible magnetic and orbital 
orderings (section^) and investigate the electronic struc- 
ture from our ab initio calculations (section IIII|) . In sec- 
tion [W] we present the results from the standard GGA 
calculations and in section we try to incorporate elec- 
tron correlation effects by means of the LDA+U approx- 
imation. 



II. STRUCTURE AND MAGNETIC 
ORDERINGS 

BaCoOs crystallizes in the hexagonal-2H pseudo- 
perovskite structure with space group P63/mmc, and lat- 
tice parameters a=5.645 A and c=4.752 A>i4 The posi- 
tion of the oxygen atoms has a free structural parameter 
whose value was taken from Ref. where it has been 
obtained from LAPW calculations. 

Figure ^ shows that BaCoOs forms chains of face- 
sharing CoOe octahedra with a much shorter (~ 2.38 A) 
Co-Co distance along the c axis than between the chains 
(~ 5.65 A), giving the structure a strong ID character. 
The Co chains are organized in a triangular manner, giv- 
ing rise to magnetic frustration when antiferromagnetism 
occurs in the plane perpendicular to the Co chains and 
assuming collinear magnetism. The Ba ion acts merely 
as an electron donor. 

In Fig. |21 we show the different magnetic orderings 
studied in this work. Besides a ferromagnetic (FM) or- 
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FIG. 1: (Color online) Hexagonal structure of BaCoOs. Ob- 
serve the chains of face-sharing C0O6 octahedra and the trian- 
gular arrangement of these chains forming a hexagonal struc- 
ture. 




FIG. 2: (Color online) Different magnetic configurations of 
BaCo03 studied. Only Co atoms are shown, dark (light) 
spheres indicate spin up (down) character, a) FM; b) A-type 
AF (FM in plane, AF between planes); c) AF-type I (AF 
along c); d) AF-type II (FM along c). 



dering (Fig. Et) also an "A-type" antiferromagnetic (AF) 
order is possible, where FM hexagonal planes are coupled 
antiferromagnetically (Fig. The triangular geome- 

try of the Co chains leads to magnetic frustration of any 
possible collinear AF long-range order in the hexagonal 
plane (perpendicular to the chains)^ Nevertheless, we 
have studied two partly AF arrangements, where each 
Co has 2 FM and 4 AF neighbors within the hexagonal 
plane. These planes can be coupled antiferromagnetically 
(AF-type I, (Fig. |2t)) or ferromagnetically (AF-type II, 
(Fig. OJi)) along the c axis. 

We study these different configurations mainly in or- 
der to explore the relative strength of inter-chain with 



respect to the intra-chain magnetic interactions. The 
quasi-one-dimensional structure of the compound would 
imply that the intra-chain magnetic coupling is most im- 
portant, since the corresponding Co-Co distance is more 
than two times smaller than the interchain Co-Co dis- 
tance. 



III. COMPUTATIONAL DETAILS 

All calculations in this paper were carried out with 
the WIEN2k package.!^ This uses the full-potential 
APW+lo methodii that makes no shape approximation 
to the potential or density. Local orbitals (Co 3s, 3p, 3d, 
Ba 4d, 5s, 5p, O 2s and 2p) were added to improve the 
flexibility of the basis set and to describe semicore states. 
The values of the atomic sphere radii (R m t) were chosen 
as: 1.70 a.u. for Co, 2.0 a.u. for Ba and 1.55 a.u. for O. 

We calculated the properties of the system within both 
the generalized gradient approximation (GGA) using the 
PBE scheme^ and the LDA+U approximation>i2i2£- In 
the latter, the DFT (orbital independent) part of the 
potential was treated by the same GGA scheme. For 
GGA, we explored the energy vs. magnetic moment in 
the range S=0 to S=l/2 using the fixed spin moment 
(FSM) methodic that fixes the total magnetic moment 
in the unit cell. We used a plane wave cutoff described by 
RmtK maa: =7 and 72 k-points in the irreducible Brillouin 
zone (IBZ) (1000 k-points in the whole Brillouin zone). 
Convergence was checked up to R m tK maa; =8 and 150 k- 
points in the IBZ. For the larger AF unit cells, equivalent 
k- meshes were employed. 

The use of the LDA+U method is especially suited for 
this moderately correlated transition metal oxide, since it 
allows to introduce an on-site Coulomb repulsion term, 
characterized by a Hubbard U. In addition, the intro- 
duction of an orbital dependent potential allows to study 
orbital orderings, which turn out to be an important issue 
in this class of compounds. 

IV. GGA CALCULATIONS 

A. FSM curve 

In order to investigate the magnetic properties of the 
present compound, we analyzed the energy of the FM 
system by varying the magnetic moment of the cell close 
to the experimentally found low-spin (LS) state. The 
unit cell contains 2 Co ions, and thus the expected total 
magnetic moment for a perfect LS state of Co would be 
2 fis- For this purpose we used the FSM method within 
the GGA approximation for the exchange-correlation po- 
tential. This method allows to calculate the total energy 
as a function of the fixed total magnetization in the unit 
cellia 

The results are shown in Fig. where the minimum 
of the FSM curve is located near (at around 1.5 \xb per 
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FSM curve (Energy vs. Magnetic Moment) 




Total magnetic moment per unit cell ( fl B ) 



FIG. 3: Energy vs magnetic moment by FSM calculation us- 
ing GGA. The minimum is near 1.5 fiB while the perfect LS 
state (2 (ib) is comparable in energy to the nonmagnetic state. 
The line is just a guide for the eye. 
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FIG. 4: Total DOS curves within GGA for fixed values of the 
total magnetic moment per unit cell of and 2/is, with S=0 
and S=0.5, respectively. 



unit cell) but definitely below the LS state. Only about 
0.5 iib reside inside each Co sphere, but the rest of the 
magnetic moment resides inside the O spheres and the 
interstitial region. The LS state would imply a half- 
metallic electronic structure, i.e. a gap in the spin-up, 
but a metallic spin-down DOS, in contrast to the exper- 
imentally observed semiconducting behavior, which can- 
not be obtained within GGA for the FM case. 

Figure 01 shows how the DOS curves change when the 
total magnetic moment is varied from to 2/xb. Half- 
metallicity is reached when the magnetic moment per Co 
site is set to an integer value. The spin-down compared 
to the spin-up DOS is approximately shifted rigidly. 

These results roughly agree with LMTO-ASA 



TABLE I: Relative energies (with respect to the FM ground 
state) of the different magnetic configurations studied with 
GGA. The value of the magnetic moment is for the Co atom 
within its atomic sphere. 





Energy (meV/Co) 


Magnetic 


Moment ((j,b) 


Ferromagnetic 







0.48 


Nonmagnetic 


5.5 




0.0 


AF-type II 


5.8 




0.11 



calculations, 8 which find the ground state to be a FM 
low-spin state. They find, however, a much higher sta- 
bilization for that state with respect to the nonmagnetic 
case, which is probably due to the atomic sphere approx- 
imation (ASA). 



B. Antiferromagnetic calculations 

The nonsymmorphic symmetry operations involving a 
translation of (0,0,1/2) along the Co chains cause degen- 
erate bands that appear near the Fermi level. To get 
an insulating behavior those degenerate bands must be 
split by breaking this symmetry, e.g. with an antiferro- 
magnetic order. Within GGA calculations, the A- type 
AF order is unstable and converges to a nonmagnetic 
solution. Among the cases studied here, the only AF 
structure that converges to a solution different from a 
nonmagnetic case is the AF-type II structure (see Fig. 
|2Jl). This ordering assumes a partial AF coupling of FM 
Co chains. The in-plane coupling is FM along the b axis, 
leading to 2 FM but 4 AF nearest neighbors. This solu- 
tion has a fairly small magnetic moment on the Co sites, 
of just 0.11/z_b per Co sphere. Although the Co chains are 
quite far appart, we obtain a much smaller moment with 
this partly AF coupling than the completely FM aligned 
case (0.48 /is per Co sphere). Apparently the inter-chain 
coupling is much stronger than one would intuitively as- 
sume according to the quasi-one-dimensional structure of 
the material. The energy of the AF-type II structure is 
higher than the FM case, which is the most stable solu- 
tion within GGA, and is even higher than a nonmagnetic 
solution (see Table HJ. We see that the energy differences 
involved are fairly small. 



C. Peierls distortion 

Another symmetry breaking, which could open a gap 
around the Fermi level, to yield the experimentally ob- 
served semiconducting state, would be a dimerization of 
the Co atoms along the chains. In order to simulate the 
Peierls distortion, we displace one of the Co atoms in the 
Co chains along the c axis towards the other Co to form 
dimers, i.e one Co is displaced from the equilibrium posi- 
tion halfway between the other type of Co atoms, which 
are at the corners of the unit cell (see Fig. QJ. For this 
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constrained Co dispalacement, the positions of the oxy- 
gen and barium atoms were relaxed by minimizing the 
forces acting on them. 

Even at very high distortions, we did not observe 
a semiconducting state, contrary to the results from 
LMTO-ASA calculations 8 (probably because these au- 
thors did not fully relax the O and Ba atoms). We take 
oxygen relaxation into account and observe a flattening 
of the bands around the Fermi level, but the hybridiza- 
tion makes the bands broad enough to overlap causing 
a metallic behavior. We observe a decrease in the total 
magnetic moment per unit cell from 1.5 [Ib in the sym- 
metric case to 1.3 (J-b in the most distorted case stud- 
ied. Since the corresponding energies are high, we can 
conclude that in this system a Peierls distortion is not 
likely to occur. A small Co distortion of 0.07 A from 
the equilibrium raises the energy by 71 meV/Co but a 
larger distortion (of 0.32 A) by 1 eV/Co. In the earlier 
bandstructure study^ it has been shown that no Fermi 
surface nesting is present. 

V. LDA+U CALCULATIONS 

In previous sections we have shown that it is not possi- 
ble to find an insulating behavior within GGA, not even 
by symmetry breaking. The LDA+U method within the 
so-called "fully localized limit"— is a scheme that approx- 
imately allows to treat the electron-electron interactions 
that occur in transition metal oxides. This method is 
known to improve over GGA in these moderately corre- 
lated transition metal oxides and to predict better values 
of the magnetic moments and gaps. 7 

The method needs two input values, the on-site 
Coulomb repulsion U and the exchange energy J. The 
value of J was taken as the average energy shift between 
spin-up and spin-down bands of the GGA spin-polarized 
calculation (approximately 0.5 eV). For this type of com- 
pounds the value of U is known to be in the range of sev- 
eral eV up to about 8 eV and was assumed to be 5 eV. 
We checked how the results depend on U and found that 
for U=2.5 eV, we basically recover the GGA solution. 
For bigger values of U (6 eV and 7 eV) , the FM solution 
remains stable with a total magnetic moment of 1.00 [Ib 
per Co ion, but the energy difference between the FM 
and AF solutions gets gradually reduced as U increases. 
Also, the magnetic moments get more localized inside the 
atomic spheres as U increases, as expected. 



A. Possible magnetic and orbital orderings 

We investigated the different magnetic configurations 
shown in Fig. [fusing the LDA+U method. In this case, 
a magnetic solution is possible for all structures, with 
strongly localized spin moments on the Co sites and neg- 
ligible contributions from the other spheres. The main 
results are summarized in Table ITT1 



TABLE II: Main results for the different magnetic configura- 
tions studied within LDA+U (U=5eV). Orbital order labelled 
"alt. c" means that the hole is alternating between A x 2 _ y 2 
and d xy orbitals in Co atoms along the Co chain, "alt. ab" 
means that it is alternating within the hexagonal plane. The 
magnetic moments M correspond to the amount of the mo- 
ment that resides inside the atomic spheres specified. 



Magnetic 


Orbital 


Energy 


M(Col) M(Co2) Insulator 


order 


order 


(meV/Co) 


(Ms) 


(Ms) 






alt. c 





0.846 


0.850 


Yes 


FM 


d a .2_ y 2 


71 


0.698 


0.698 


No 




d x y 


76 


0.705 


0.705 


No 




alt. c 


7 


0.688 


-0.689 


Yes 


A-type AF 


d x 2_ y 2 


52 


0.768 


-0.767 


Yes 




d X y 


54 


0.770 


-0.770 


Yes 


nonmagnetic 


none 


58 








No 




alt. c 


12 


0.811 


-0.818 


Yes 


AF-type I 


d x 2_ y 2 


58 


0.897 


-0.897 


Yes 




d x y 


57 


0.900 


-0.900 


Yes 




alt. ab 


86 


0.651 


-0.674 


No 


AF-type II 




93 


0.683 


-0.683 


No 




d x y 


93 


0.682 


-0.682 


No 




FIG. 5: (Color online) Spin density plot (isosurface at 0.2 
e/A 3 ) for the alternating orbital ordered A-type AF configu- 
ration. The spin moment is located in A x 2_ y 2 and d xy orbitals 
alternating along the Co chains, (the orbitals form angles of 
45 degrees with each other). The FM orbital-ordered solu- 
tion would be equivalent but with all Co in equal color (same 
spin). 



Co 4+ in a (nearly) octahedral environment would be a 
t|g e g system with one hole in the spin-down t2 g states. 
In the present hexagonal crystal structure, these are the 
d 2 2, d xy and d x 2_ y 2 orbitals, whereas the d xz and d yz 
orbitals form the e g manifold pointing towards the O lig- 
ands. It turns out that within LDA+U the hole can be 
localized predominatly in the d x 2_ y 2 and/or d xy orbitals, 
but never in the d z 2 state, which is always fully occupied. 
From Table m we see that the energy difference between 
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solutions with a d x 2_ y 2 or d^j, hole is very small, but both 
are much higher than the solution with an alternating 
orbital ordering (i.e. Co atoms along the Co chain have 
alternating d x 2_ y 2 and d xy holes). This energy gain is 
much larger than the differences between different mag- 
netic order (FM vs. A- type AF or AF-type I or even 
a nonmagnetic solution). Therefore, magnetic configu- 
rations play a minor role in this quasi-one-dimensional 
system. This finding is consistent with the results in 
Ref. |^] where it was shown that this type of orbital 
ordering is caused by superexchange interactions in a 
frustrated Jahn- Teller system with metal-oxygen-metal 
bonds of 90°. In BaCoC>3, this angle is 86.9 and thus 
spin and orbital degrees of freedom are decoupled, mak- 
ing the spin exchange much weaker than the orbital one. 

The A-type AF structure is the most stable antiferro- 
magnetic configuration studied here, ft is the only one 
that can lead to a collinear long-range AF ordering (along 
the Co chains). As mentioned above, the energy differ- 
ence to the FM ground state decreases with increasing 
U. 

In the AF-type II case we restricted the Co atoms to 
be equivalent within the chain, and allowed orbital order 
only between the chains. Even this case of inter-chain 
orbital ordering is energetically favourable (though not 
as much as within the Co chain). Once again this is a 
signal for sizable inter-chain interactions, which are not 
expected to be so large in a quasi-one-dimensional sys- 
tem. 

For the FM case, the introduction of an alternating 
orbital ordering along the chains removes the symmetry 
that makes the two types of Co atoms along c equivalent 
and thus allows to produce an energy gap around the 
Fermi level, leading to an insulating behavior. Such an 
insulating state can also be reached when this symmetry 
is broken by an AF ordering along the chain. In the case 
of AF-type II, the intra-chain symmetry is not broken 
and thus a metallic state is found. Just as in GGA, the 
FM solution is the most stable in LDA+U calculations. 
In the latter case, the total magnetic moment within the 
unit cell (which now contains 4 Co atoms) is 4.00 /is 
for the alternating orbital ordering (which corresponds 
to the LS state), while it is only 3.9 hb for the other 
orbital orderings (with metallic character). 

Figure [S] shows the spin density of the alternating or- 
bital ordered state of the A-type AF structure. The FM 
case would look similar, but the colors (indicating spin- 
up and spin-down) would be equal for every Co atom. 
The spin-density along the Co chains is localized in two 
orbitals, which are rotated 45 degrees to each other, hav- 
ing their lobes in the hexagonal plane (d x 2_ y 2 and d xy , 
respectively). In addition we see some spin density on the 
O atoms indicating the importance of the 90° superex- 
change mentioned above. 
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FIG. 6: Total density of states (DOS) for the different mag- 
netic configurations in the alternating orbital ordered state. 
Semiconducting behavior is related to symmetry breaking 
along the Co chains. The energy is taken with respect to 



B. Densities of states (DOS) 

In this section we discuss the main features of the elec- 
tronic structure of the system using the DOS plots. In 
GGA, the total moment of 2 {1b per unit cell is fairly 
delocalized over the entire unit cell and thus leads to a 
relatively small moment inside the Co spheres (see Ta- 
ble HJ). Within the LDA+U method, we obtain a phys- 
ically more plausible solution, since magnetic moments 
are more localized in the Co spheres, as expected for a 
correlated system (see Table [nj. Also, we get values of 
the magnetic moment in agreement with the experimen- 
tal LS state and the "hole" (see above) is localized in one 
of the t2 g -like orbitals. Using the GGA approximation we 
cannot obtain the experimentally found semiconducting 
state, not even when we break the symmetry by a Peierls 
distortion. This semiconducting state, however, can be 
obtained by means of the LDA+U approach. 

In Fig. we show the total DOS for the four mag- 
netic configurations studied. All of them are in an al- 
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FIG. 7: Partial Co d DOS for the Col atom in a ferromag- 
netic, alternating orbital ordered case. The hole is located at 
the d x 2_ y 2 orbital. 



FIG. 8: Partial Co d DOS for the Co2 atom in a ferromag- 
netic, alternating orbital ordered case. The hole is located at 
the d xy orbital. 



ternating orbital ordered configuration, since this is their 
lowest energy state. The main observation is that insu- 
lating behavior is found for the FM, the A-type AF and 
AF-type I structures, where alternating orbital ordering 
occurs along the Co chains. In the AF-type II case we 
forced the translational symmetry along the Co chains 
and thus excluded intra-chain orbital ordering (we have 
FM chains coupled AF to some of their neighbors). In 
this case, no splitting appears in the bands around the 
Fermi level. Hence, the lowest energy states with an al- 
ternating orbital ordering reproduce the experimental ev- 
idences, while without this symmetry breaking a metallic 
solution would occur. 

In the case of FM coupling and alternating orbital or- 
dering, we show in Figs. and [S] the partial DOS corre- 
sponding to the different Co d orbitals for the two types 
of adjacent Co atoms along the Co chains. For Col, the 
hole is mainly located in the d x 2_ y 2 orbital with some 
d xy and d xz component (see peaks just above Ef), where 
the latter produces the tilting visible in the spin-density 
shown in Fig. |SJ In the case of Co2, the hole is mainly 
of d xy character with some admixture of d yz and d x 2_ y 2. 



The e g -like orbitals (d^ and d yz ) are predominantly un- 
occupied while the d z 2 orbital is fully occupied in any 
case. 

In Fig. [5] we observe that the hole found in the total 
DOS plot has some O p character as well. This is again a 
manifestation of the importance of the 90° superexchange 
mechanism. It is mainly due to the p^ and Pj, but not 
due to the p z orbital which points along the chain axis. 
One can see this from Fig. El as well, noticing that the O 
spin cloud is mainly in the plane perpendicular but not 
along the Co chains. 



VI. SUMMARY 

In this paper we carried out ab initio density functional 
theory calculations using the full-potential APW+lo 
method as implemented in the WIEN2k code on the tran- 
sition metal oxide BaCo03. 

The calculations within GGA to the exchange- 
correlation potential gave a ferromagnetic configuration 
as the most stable one. The value of the spin magnetic 
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tal observations. 

Our electronic structure calculations indicate that the 
inter-chain couplings are weak but stronger than ex- 
pected for a nearly ID system. 
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FIG. 9: Partial O p DOS for the ferromagnetic, orbital or- 
dered case. Observe that the hole has some O p character, 
which comes from p x and p y orbitals but has no component 
along the Co chain axis. 



moment was calculated using the FSM method and the 
result was 1.5 yUs/unit cell, smaller than the expected 
2 yu^/unit cell for a perfect LS state. Antiferromagnetic 
configurations within the Co chain could not be obtained. 
All the GGA solutions are metallic in contrast to the ex- 
perimental semiconducting state. We investigated the 
possibility of a Peierls distortion as a cause for the semi- 
conducting behavior, but had to rule it out by total en- 
ergy calculations. 

Within the LDA+U method, we could incorporate the 
correlation effects among the Co 3d electrons. We stud- 
ied different magnetic and orbital orderings and found 
that an intra-chain alternating orbital ordering with FM 
coupling is the most stable solution. The orbital order- 
ing is caused by a nearly 90° superexchange mechanism 
of the face-sharing octahedra and leads to the main sta- 
bilization, whereas the magnetic order is of secondary 
importance. In addition, this orbital ordering produces 
an insulating solution, which reproduces the experimen- 
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